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ABSTRACT 

We present new photometry of HD 149026 spanning five transits of its "super- 
Neptune" planet. In combination with previous data, we improve upon the determina- 
tion of the planet-to-star radius ratio: Rp/R^ = 0.0491]']qqqq5. We find the planetary 
radius to be 0.71 it 0.05 -Rjup, in accordance with previous theoretical models invoking 
a high metal abundance for the planet. The limiting error is the uncertainty in the 
stellar radius. Although we find agreement among four different ways of estimating 
the stellar radius, the uncertainty remains at 7%. We also present a refined tran- 
sit ephemeris and a constraint on the orbital eccentricity and argument of pericenter, 
ecosw = —0.0014 it 0.0012, based on the measured interval between primary and sec- 
ondary transits. 

Subject headings: planetary systems — stars: individual (HD 149026) 

1. Introduction 

Many clues about the processes of planet formation and evolution have been discovered by 
studying the ensemble properties of exoplanets, such as the "brown dwarf desert" (Halbwachs et 
al. 2000, Marcy & Butler 2000) and the tendency for metal-rich stars to have more detectable 
planets (Santos et al. 2003, Fischer & Valenti 2005). However, there are also individual exoplanets 
whose properties bear directly on theories of planet formation and evolution. One of the best 
examples is the transiting planet HD 149026b (Sato et al. 2005). 

Compared to Saturn, HD 149026b has a similar mass but its radius is 15% smaller, despite the 
intense irradiation from its parent star that should enlarge the radius. Sato et al. (2005) modeled 
HD 149026b as a dense heavy-element core surrounded by a fluid envelope of solar composition. 
They found a core mass of 70-80 M^, which is 65-75% of the total mass of the planet. This is larger 
than the canonical core mass of 10-20 that is expected from the core-accretion theory of planet 
formation (Mizuno 1980, Pollack et al. 1996). The finding of a highly metal-enriched composition 
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was confirmed in models by Fortney et al. (2006), Ikoma et al. (2006), Broeg &: Wuchterl (2007), 

and Burrows et al. (2007). The latter authors dubbed HD 149026b a "super-Neptune" because the 
inferred mass fraction of heavy elements is similar to that of an ice giant rather than a gas giant. 

Interestingly, the parent star has a rather high metallicity ([Fe/H] = +0.36; Sato et al. 2005). 
The observation of a large core in such a metal-rich system would seem to support the core-accretion 
theory as opposed to coreless alternatives such as gravitational instability (Boss 1997). However, 
the larger-than-expected core mass raises some questions. Why did the growing protoplanet not 
accrete gas efficiently? Or if it did, what happened to its envelope of light elements? Many scenarios 
have been proposed: a collision of two massive protoplanets (Sato et al. 2005, Ikoma et al. 2006), 
in situ formation in a low-pressure nebula (Broeg & Wuchterl 2007), a viscous and evaporating 
gas disk (Ikoma et al. 2006), and a separation of gas from planetesimals at the magnetospheric "X 
point" (Sato et al. 2005). 

More recently, Harrington et al. (2007) found that the 8 fim brightness temperature of HD 149026b 
exceeds its expected blackbody temperature, even if the planet is assumed to absorb all of the inci- 
dent stellar radiation. In this sense the planet is anomalously hot. The high temperature may result 
from novel atmospheric or structural properties. Most recently, Torres et al. (2007) announced the 
discovery of a transiting planet, HAT-P-3b, whose measured mass and radius indicate that it too 
is highly enriched in heavy elements. 

In short, HD 149026b seems to be the harbinger of an entirely new kind of planet that current 
models of planet formation, evolution, and structure cannot accommodate without interesting and 
possibly exotic modifications. Because of this situation, it is desirable to improve the reliability 
and the precision of estimates of the system parameters, and especially a key parameter that makes 
this planet unusual: its small radius. 

One can measure the planetary radius by gathering photometry during transits, modeling 
the light curve, and supplementing the model with external information about the stellar radius. 
Previously, Sato et al. (2005) analyzed 3 light curves, and Charbonneau et al. (2006) added 3 light 
curves. In this paper we present another 5 light curves of comparable or higher quality to the 
previously published data, and we simultaneously model all of the data to derive the most precise 
planetary, stellar, and orbital parameters that are currently available. We present our observations 
and data reduction procedure in § 2 and the light-curve modeling procedure in § 3. We provide 
the results in § 4, along with an extended discussion about the limiting error: the uncertainty in 
the stellar radius. The final section summarizes the results and speculates on future prospects for 
improvement. 



2. Observations and Data Reduction 



We used three of the 0.8 m automated photometric telescopes (APTs) at Fairborn Observatory 
to measure the transits of HD 149026b that occurred on UT 2006 April 26, 2006 May 20, 2007 May 3, 
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2007 June 18, and 2007 June 21. We observed the first three transits with the Til APT and 
observed the last two transits simultaneously with the T8, TIO, and Til APTs. All three telescopes 
are equipped with two temperature-stabilized EMI 9124QB photomultiplier tubes for measuring 
photon count rates simultaneously through Stromgren b and y filters. 

On a given night, each telescope automatically acquired brightness measurements of HD 149026 
{V = 8.15, B -V = 0.61) and the comparison star HD 149504 {V = 6.59, B -V = 0.44), which 
was previously demonstrated to be stable in brightness at the 0.002 mag level or better (Sato et 
al. 2005). We also measured the dark count rate and the sky brightness in the vicinity of each star. 
We used a diaphragm of diameter 45" for all the integrations. The integration time was 20 seconds 
on the comparison star and 30 seconds on the (fainter) target star. We computed the magnitude 
difference for each pair of target-comparison observations. To increase the signal-to- noise ratio of 
each measurement, the differential magnitudes from the b and y pass bands were averaged, resulting 
in a differential magnitude for a synthetic {b+y)/2 pass band. The typical cadence of the differential 
magnitude measurements was 1.4 min for the first 3 transits and 1.2 min for the last 2 transits. 

For each raw light curve, we fitted a linear function of time to the out-of-transit data and 
divided the data by this function. This was intended to correct for differential airmass effects and 
other systematic errors in the photometry, and also to normalize each light curve to have unit 
mean flux outside of the transit. (We also tried fitting a function of airmass rather than time, 
but this gave slightly poorer results.) The final light curves are shown in Fig. 1, and the data 
are given in Table 1. The standard deviation of the out-of-transit data is approximately 0.2% in 
all cases, which is typical for APT observations of bright stars. For additional information on the 
telescopes, photometers, observing procedures, data reduction techniques, and typical photometric 
precision, see Henry (1999) or Eaton, Henry, & Fekel (2003). The lower right panel in Fig. 1 is a 
composite light curve of all 8 transits observed by the APTs from Sato et al. (2005) and this paper. 
The composite light curve was created by subtracting the mid-transit time from each of the time 
stamps, and then averaging into 30 second bins. It is shown here for display purposes only; the 
fitting procedure described in the next section was carried out on the unbinned data. 



3. Determination of System Parameters 

We fitted all of our new photometric data jointly with the three [h+y) /2 light curvespresented 
by Sato et al. (2005) and the g and r light curves presented by Charbonneau et al. (2006) o We used 
a parameterized model based on a two-body circular orbit. The orbit is specified by the masses 
of the star and planet {M^ and Mp), the inclination with respect to the sky plane (z), the orbital 
period (P) and a particular midtransit time (Tc). The star and planet are taken to be spheres 



^We did not include the l/-band light curve of Charbonneau et al. (2006) because of its comparatively large errors 
and sparse time sampling. 
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Fig. 1. — Stromgren (6 + y)/2 photometry of five new transits of HD 149026, acquired with three 
of the 0.8 m APTs at Fairborn Observatory. The lower right panel is a composite hght curve 
created from the five hght curves reported here as well as the three {b + y)/2 light curves previously 
published by Sato et al. (2005). 
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with radii i?* and Rp respectively, and when their sky-projected centers are within Ri, + Rp oi 
one another we use the Mandel & Agol (2002) formulas to compute the flux decrement due to the 
partial blockage of the limb-darkened stellar surface. This is the same code that has been developed 
for the Transit Light Curve project (Holman et al. 2006, Winn et al. 2007). For HD 149026, we 
assumed the limb darkening law to be linear, with a coefficient given by Claret (2000) for a star of 
the appropriate temperature and surface gravity. For the {b + y)/2 data we used a coefficient of 
0.712, which is the mean of the tabulated b and y coefHcients. 

Not all of the parameters listed above can be determined from transit photometry alone. One 
set of parameters that can be determined from an individual light curve is Tc, Rp/Ri,, a/R^, and i, 
where a is the semimajor axis. Our approach was to fix Af^, Adp, and P at previously determined 

values (thereby fixing a through Kepler's third law), and then fit for Rp, i?^, i, and Tc- The results 

1/3 1/3 

for R^ and Rp are specific to the choice of M^, but they scale as because a oc M^' when the 

uncertainty in P is negligible, as it is here. We assumed = 1.3 Mq, following Sato et al. (2005), 
a choice that was subsequently corroborated by our analysis of the observable stellar properties 
and the results for a/R^ (see § 4.2). 

The fitting statistic was 



where Nf is the number of flux measurements, /j (obs) is the flux observed at time j, aj controls 
the weights of the data points, and /j(calc) is the calculated flux. Experience has shown that the 
data weights aj should account not only for the single-measurement precision but also the time- 
correlated ("red") noise that afflicts most time-series photometry (see, e.g., Gillon et al. 2006). 
The most important timescale in a transit light curve is the ~10 min duration of the ingress and 
egress, since the resolution of ingress and egress is what permits the determination of a/Ri, and 
i in addition to Rp/R^. To assess the noise on this timescale, we first calculated the standard 
deviation of the unbinned out-of-transit data ((Ji) for each light curve. Then we averaged the out- 
of-transit data into 10 min bins consisting of N data points, where N depended on the observing 
cadence, and recalculated the standard deviation ((Tat). In the absence of red noise, one would 
expect (Ttv = (Ti/^/N, but in practice fijv was larger than ai/^/N by some factor /3. Therefore, we 
set the data weights equal to (3 ai. The results for f3 ranged from 1.05 to 1.27. 

We used a Markov Chain Monte Carlo algorithm to determine the best-fitting parameter values 
and confidence intervals. This algorithm delivers an estimate of the a posteriori joint probability 
distribution for all of the parameters (see Holman et al. 2006 or Winn et al. 2007 for more details). 
For each parameter, we took the mode of the distribution after marginalizing over all other pa- 
rameters to be the "best value." We defined the 68% confidence limits pio and phi as the values 
between which the integrated probability is 68%, and for which the two integrals from pmin — P\o 
and Phi Pmax were equal. 




(1) 



At first, the preceding computations were performed using free parameters for both the orbital 
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period and a single midtransit time. Thus the transits were required to be spaced by integral 
multiples of a fixed period. However, we also wanted to measure the individual transit times, in 
order to search for variations that might be indicative of additional bodies in the planetary system 
(Agol et al. 2005, Holman & Murray 2005). To do this, we fixed Rp, i?*, and i at the best values 
determined in the first step, and then performed a three-parameter fit of each individual light curve. 
The parameters were Tc along with the zero point and slope of the linear function that was used to 
correct the out-of-transit data. (Fixing the values of Rp, i?*, and i is justified because the errors 
in those parameters are not correlated with the error in the transit time.) We did this not only for 
the 5 new light curves, but also for the 5 previously published light curves, to provide consistency 
in the treatment of errors. We then used these transit-time measurements to refine the estimates 
of P and Tc (see § 4.3). For our final results for the photometric parameters Rp/R^,, a/Ri,, and 
i, we reran the MCMC algorithm on the entire data set, using fixed values of P and Tc from our 
refined ephemeris. 



4. Results 

The results of the light-curve analysis are given in Table 2 and discussed in § 4.1. The transit 
times of the individual light curves are given in Table 3, and in § 4.2 we use those times to derive a 
new transit ephemeris. We also use the new ephemeris along with a previously-measured midpoint 
of a secondary eclipse to place an upper bound on one aspect of the orbital eccentricity. In order 
to derive the actual planetary radius (as opposed to the planet-to-star radius ratio), one must 
supplement the light-curve analysis with external information about the stellar radius or mass. In 
§ 4.3 we investigate four different methods to estimate the stellar radius, the results of which are 
given in Table 4. The planetary parameters derived from the light curves and the stellar radius 
estimates are given in Table 5. 



4.1. Photometric Parameters 

The a posteriori probability distributions for Rp/R^,, a/Ri,, and i are shown in Fig. 2. The 
most well-constrained of these three basic light curve parameters is the radius ratio, Rp/R^ = 
0.04911q qqq5, with a precision of approximately 2.3%. The radius ratio is determined largely by 
the observed transit depth, which is the smallest among all of the 20 transiting planets known to 
date. While there are smaller planets, such as the Neptune-sized GJ 436b (Gillon et al. 2007), they 
all orbit smaller stars, making their radius ratios and transit depths larger than that of HD 149026. 

Less well-constrained are a/i?^ and i, the parameters that depend on the observed durations 
of the ingress, egress, and total phases of the transit. Table 2 gives the results for those parameters 
as well as the impact parameter and the durations, which can be derived in terms of Rp/Ri,, 
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Fig. 2. — Estimated a posteriori probability distributions from the joint fit to the transit hght 

curves. The top panels show the single- variable distributions, in which the mode is marked with 
a solid line and the 68% confidence limits with dashed lines. The bottom panels show the two- 
dimensional distributions, in which the contours mark the 68% and 95% confidence limits. 
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and i. The data are consistent with impact parameters ranging from —0.36 to 0.36. As with 

any echpsing binary system, the data cannot distinguish between positive and negative impact 
parameters; the probability distributions are perfectly symmetric about b = and i = 90°. In this 
case the peak probability occurs at the central values b = 0, i = 90°. The quantity a/R^ has a 
highly asymmetric error bar. The observed duration of the entire transit event enforces the upper 
limit on a/Ri,, while the ratio of the ingress (or egress) duration to the total duration enforces the 
lower limit on a/Ri,. 

Our findings are consistent with the two previous light-curve analyses, by Sato et al. (2005) 
and Charbonneau et al. (2006), although our analysis method is different in several ways besides 
the use of an expanded dataset. First, we have attempted to account for time-correlated noise in 
the photometry, which was neglected in the previous analyses. Second, unlike the previous authors, 
we have not incorporated any a priori constraints on the stellar properties into our fitting statistic. 
We made this choice in order to clarify what information is derived from the light curves themselves; 
for example, the previous works did not call attention to the results for Rp/R^ even though that 
parameter is more precisely known than either Rp oi R^. In addition, our analysis method provides 
an estimate of a/R^, that is independent of any assumptions about the parent star, except for the 
very weak dependence on the chosen limb darkening parameter. This is useful because a/-R^ can 
be used to determine the stellar mean density (Seager & Mallen-Ornelas 2003, Sozzetti et al. 2007, 
Holman et al. 2007): 

The last term in this expression may be neglected in this case because pp ~ pi, and {Rp/Ri,)^ ~ 10~^. 
Our independent estimate of is useful in characterizing the parent star, as described below. 



4.2. The Stellar radius 

To determine the quantity of intrinsic interest, Rp, we can multiply our result for Rp/R^, by a 
value of Ri, obtained by other means. We have investigated four diflFerent methods for determining 
R*: 

Stefan-Boltzmann Law. — The bolometric luminosity, effective temperature, and photospheric 
radius of HD 149026 are related via Lboi = 47ri?^(TrJj. We use the Hipparcos parallax and apparent 
magnitude (tt = 12.68 ± 0.79 mas, V = 8.15 ± 0.02; Ferryman et al. 1997) to compute the absolute 
V magnitude, apply a bolometric correction of —0.027 it 0.014 (Flower et al. 1996), and use the 
spectroscopically determined T^s = 6147±50 K (Sato et al. 2005). The result isR^ = 1.46±0.10 Rq. 
This is essentially identical to the value quoted by Sato et al. (2005) who used the same method. 

Spectral Energy Distribution Fit. — Masana et al. (2006) presented an alternative means of 
estimating the effective temperature and bolometric correction, using VJHK photometry. They 
also provided radius estimates for many nearby stars based on this technique. Using the Hipparcos 
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parallax and V magnitude along with 2MASS near-infrared photometry, their result for HD 149026 
is = 1.515 ±0.096 Rq. 

Yonsei-Yale Isochrone Fit. — Stellar evolutionary models may be used to estimate the mass, 
radius, and age of a star with a given effective temperature, luminosity (or gravity), and metallicity. 
We used the Yonsei-Yale models (Yi et al. 2001, Demarque et al. 2004) because they are conveniently 
provided with tools for interpolating isochrones in both age and metallicity. For the effective 
temperature, we used T^s = 6160 ± 50 K, a weighted mean of the results of Sato et al. (2005) and 
Masana et al. (2006). We used the photometric result for a/i?^ as our our proxy for surface gravity, 
and we explored the rang^of metallicities [Fe/H] = 0.36 it 0.08. For each metallicity we considered 
a range of ages from 0.1 to 14 Gyr, in steps of 0.1 Gyr. We interpolated the isochrones using a fine 
mass grid and compared the points with the measured values of T^s and a/Ri,. We computed 
at each point based on the modeled and observed values of Teff, a/R^,, and metallicity. Then we 
weighted the points by exp(— x^/2) and applied an additional weighting to take into account the 
density of stars on each isochrone, assuming a Salpeter initial mass function. The "best-fitting" 
stellar properties were taken to be the weighted mean of the properties of all the points. For more 
details and other applications of this analysis, see Torres, Winn, & Holman (in preparation). For 
HD 149026, the results are = 1.294to:o5o Mq, = 2.430;[]|^| L©, and R^ = IMStoUt ^Q- 
Similar results were obtained when the spectroscopically determined value of log g was used instead 
oi a/R^. The theoretical isochrones and the observational constraints are shown in Fig. 3. Kepler's 
Law with Stellar Mass Prior. — As mentioned earlier, the quantity a/i?^ that is determined from 
the transit photometry can be used to find (Eq. 2). With an a priori estimate of M^, one may 
use to determine i?^,. Taking M^, = 1.30 it 0.06 Mq based on the isochrone fit described above, 
we find R^ = l-SbtoH 

All of the results for the stellar radius are summarized in Table 4. They are all consistent with 
one another at the la level, with a weighted mean of 1.45 Rq. However, it must be emphasized 
that while the methods are different, they are not wholly independent. The first two methods both 
rely on the Hipparcos parallax, which is the largest source of error in both cases. The latter two 
methods both rely on the Yonsei-Yale stellar evolutionary models. For this reason we cannot say 
confidently that the uncertainty in i?^ is any smaller than the uncertainty in each of the individual 
measurements, although the mutual agreement is certainly reassuring. In what follows we adopt 
the consensus value i?^, = 1.45 it 0.10 Rq, the same value used in the previous light curve analyses. 

Assuming a Gaussian error distribution for i?^, and the error distribution for Rp/ R^ obtained 
from our light curve analysis, we find the planetary radius to be Rp = 0.71 it 0.05 -Rjup. This 
can be compared to the previously published results of 0.725 it 0.050 i?jup (Sato et al. 2005) and 
0.726 lb 0.064 i?jup (Charbonneau et al. 2006), keeping in mind our different method of analysis 



^ Sato et al. (2005) reported a metallicity of [Fe/H] — +0.36 with an internal uncertainty of 0.05. To be con- 
servative, we adopted a somewhat larger uncertainty of 0.08, recognizing that different methods for determining the 
metallicity often produce systematic differences of this size. 
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Fig. 3. — Model isochrones from the Yonsei-Yale series by Yi et al. (2001) and Demarque et 
al. (2004), corresponding to ages of 1-14 Gyr (left to right), for the measured composition of 
[Fe/H] = +0.36, along with the observational constraints. Top: The vertical axis is log 5, and the 
shaded box shows the la range based on the spectroscopically determined value of log g. Bottom: 
The vertical axis is a/R^, which is proportional to the cube root of the stellar mean density (see 
Eq. 2). The shaded box shows the la range based on the photometrically determined value of 
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and treatment of observational errors. The results are all in agreement. Indeed the differences are 

smaller than one would expect from Gaussian statistics, given the quoted error bars, though wc 
note that 5 of the 10 light curves we fitted were taken from those previous works. The precision in 
Rp is not improved because the limiting error is the uncertainty in R^, which is unchanged. 



4.3. Transit Times 

For planning future observations of this system it is important to be able to predict transit 
times as precisely as possible. We used all of the transit times given in Table 3 to calculate a 
photometric ephemeris for this system, 

n(E) = n{o) + EP, (3) 

where Tc is the transit midpoint, E is the integral transit epoch, and P is the orbital period. The 
linear fit had x^/^doi = 0.63 and A^dof = 9) suggesting that the errors quoted in Table 3 have been 
somewhat overestimated. The results are: 

Tc(0) = 2454272.7301 ± 0.0013 [HJD] (4) 
P = 2.8758882 ± 0.0000061 days. (5) 

Our value for the orbital period is in agreement with the previously published values and is about 
25 times more precise. Figure 4 is the O— C (observed minus calculated) diagram for the transit 
times. 

For a circular orbit, successive transits and secondary eclipses should be spaced by exactly half 
an orbital period. Recently, Harrington et al. (2007) observed a secondary eclipse of HD 149026 
with the Spitzer Space Telescope, allowing the assumption of a circular orbit to be checked. In the 
presence of a small but nonzero orbital eccentricity, the time difference between the midpoint of 
secondary eclipse, Tgec, and the time of transit, Ttra, is 

Tsec - Ttra ~ ^ ^1 + ^eCOSW^ , (6) 

where u is the argument of pericenter (Kallrath Sz Milone 1999, p. 62). Harrington et al. (2007) 
measured the midpoint of a secondary eclipse to be HJD 2453606.960 ± 0.001, represented by the 
open circle in Figure 4. The timing offset of equation (6) is —3.6 ±3.1 minutes, corresponding 
to ecosu = —0.0014 ± 0.0012. The orbit does indeed appear to be nearly circular, as one would 
expect from the dissipative effects of stellar and planetary tidal interactions. 



5. Discussion and Summary 



We have presented 5 new transit light curves of the exoplanet HD 149026b and analyzed them 
along with 5 previously published light curves. The joint analysis has resulted in much more precise 
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Fig. 4. — Transit timing residuals for HD 149026b. The calculated times, using the ephemeris 
derived in § 4.3, have been subtracted from the observed times. The filled symbols represent 
observations of transits. The open symbol represents the observation of the secondary eclipse by 
Harrington et al. (2007). The secondary eclipse datum was not used in the fit. 
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determinations of the orbital period and transit ephemerides, and also in a more precise value of 

the planet-to-star radius ratio. In some cases, this ratio is of primary interest, such as inferring the 
brightness temperature of the planet from the depth of a secondary eclipse (Harrington ct al. 2007), 
or testing for any wavelength dependence in the radius ratio as a means of identifying planetary 
atmospheric features (see, e.g., Charbonneau et al. 2002). 

However, when it comes to understanding the interior structure of the planet, the quantity 
of primary interest is Rp itself, and here we can offer no significant improvement. The limiting 
error is the 7% uncertainty in the stellar radius. This error was not reduced by acquiring more 
light curves, although we did find agreement between the results of 4 different (and intertwined) 
methods for estimating the stellar radius using all of the available data. Thus, we leave unchanged 
the interpretation of this planet as a being unexpectedly small for its mass, and likely to be highly 
enriched in heavy elements (Sato et al. 2005, Fortney et al. 2006, Ikoma et al. 2006, Burrows et 
al. 2006). 

Further improvement will depend upon progress in measuring the stellar radius. Baines et 
al. (2007) recently used optical interferometry to measure the angular diameter of the planet-hosting 
star HD 189733, and combined it with the Hipparcos parallax to measure the stellar radius. For 
HD 149026, similar observations are not likely to result in a more precise value of the stellar radius, 
at least not in the near future. This is not only because of the 6% uncertainty in the parallax, but 
also because the expected angular diameter is only ~180 /xas, which is only 7-8 times larger than 
the measurement error that was achieved for HD 189733. 

Supposing the parallax were known with 10 yuas precision (as one might hope from a space- 
based interferometric mission) , the error in the Stefan-Boltzmann method for determining i?^ would 
be reduced to 2.7%. The limiting errors in that case would arise from the effective temperature and 
bolometric correction. In the nearer term, a possible path forward is the continued acquistion of 
high-quality transit photometry, in order to improve upon our measurement oi a/R), and thereby 

1/3 

establish the stellar mean density with greater precision. At fixed mean density, i?^, varies as , 
and our application of the Yonsei-Yale models to HD 149026 suggests that the stellar mass has 
already been pinned down to within 4.6%. If a/i?* were known exactly, the fractional error in 
the stellar radius would be approximately 1.5% (i.e., one-third as large as the fractional error in 
the stellar mass). In effect, transit photometry measures M^/R^, and the stellar models generally 
constrain a different combination of and R^ (see, e.g., Cody &; Sasselov 2002). We encourage 
observers to be persistent in gathering additional seasons of ground-based photometry and look 
forward to the results of space-based photometry for this system. 

We are grateful to the anonymous referee for a thorough and helpful review of the manuscript. 
G.W.H. acknowledges support from NSF grant HRD-9706268 and NASA grant NNX06AC14G. GT 
acknowledges partial support for this work from NASA grant NNG04LG89G. 
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Table 1. Photometry of HD 149026 



Heliocentric Julian Date 


Relative flux 


Uncertainty 


2453852.75015 


1.0031 


0.0027 


2453852.75115 


1.0061 


0.0027 


2453852.75215 


1.0023 


0.0027 



Note. — The data were obtained with the automatic 
photometric telescopes (APTs) at Fairborn Observatory. 
Differential magnitudes were measured in the Stromgren 
b and y passbands, and the b and y results were averaged. 
The time stamps represent the Heliocentric Julian Date at 
the time of mid-exposure. The uncertainties include the 
"red noise" correction described in § 3. We intend for this 
Table to appear in entirety in the electronic version of the 
journal. An excerpt is shown here to illustrate its format. 
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Table 2. HD 149026: Transit Light Curve Parameters 



Parameter 


Value 


68% Upper Limit 


68% Lower Limit 




0.0491 


+0.0018 


-0.0005 


a/Ri, 


7.11 


+0.03 


-0.81 


i [deg] 


90.0 


+3.1 


-3.1 


h = acosi/R^ 


0.00 


+0.36 


-0.36 


Total transit duration^ 


3.254 


+0.057 


-0.028 


Ingress or egress duration*^ 0.153 


+0.052 


-0.002 



'^Defined as the time between first and fourth contacts (i.e., between the 
moments when the projected planetary and stehar disks are externally tangent). 

'^Defined as the time between first and second contacts (i.e., the duration 
over which the projected planetary disk crosses the stellar limb, from external 
tangency to internal tangency). In our model, the ingress and egress durations 
are equal. 

Note. — Results of fitting ten hght curves: three light curves [{b+y)/2] from 
Sato et al. (2005); two light curves {g and r) from Charbonneau et al. (2006); 
and five light curves [(b + y) /2] from this work. Not all of the parameters are 
independent. One may regard Rp/R^,, a/Ri,, and i as the basic parameters 
from which the other results in this table may be derived. 
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Table 3. 


HD 149026: 


Midtransit times 




Telescope 


Epoch 


Mid-transit time 


Uncertainty 




E 


[II.ID] 


[davs] 



Til 0.8m APT 


-267 


2453504. 


,8707 


0, 


,0022 


Til 0.8m APT 


-259 


2453527. 


,8732 


0, 


,0021 


FLWO 1.2m 


-259 


2453527. 


,8722 


0, 


,0022 


T8, TIO, Til 0.8m APTs 


-258 


2453530. 


,7517 


0, 


,0031 


FLWO 1.2m 


-250 


2453553. 


,7583 


0, 


,0013 


FLWO 1.2m 


-250 


2453553. 


,7598 


0, 


,0045 


Til 0.8m APT 


-146 


2453852. 


,8514 


0, 


,0024 


Til 0.8m APT 


-138 


2453875. 


,8600 


0, 


,0050 


Til 0.8m APT 


-17 


2454223. 


,8355 


0, 


,0070 


T8, TIO, Til 0.8m APTs 


-1 


2454269. 


,8554 


0, 


,0021 


T8, TIO, Til 0.8m APTs 





2454272. 


,7293 


0, 


,0017 



Note. — Based on these measurements, we derived a transit ephemeris 
Tc(^) = rc(0) + EP with rc(0) = 2454272.7301(13) [HJD] and P = 
2.8758882(61) days, where the numbers in parentheses indicate the Icr 
uncertainty in the final two digits. 



Table 4. HD 149026: Stellar Radius 



Radius Method Reference 



1.46 ±0.10 


Stefan-Boltzmann Law 


1,3 


1.515 ±0.096 


Spectral Energy Distribution Fit 


2 


1 '?fi8+0-124 
±.ouo_Q 


Yonsei-Yale Isochrone Fit 


3 


1 oc:+0.17 

±.Od_Q Q2 


Kepler's Law with Stellar Mass Prior^ 


3 



^Using M* = 1.30 ± 0.06, based on the Yonsei-Yale isochrone fit to 
a/Ri: and T^s. 



Note. — References: (1) Sato et al. 2005; (2) Masana et al. 2006; 
(3) This work. 



Table 5. HD 149026: Planetary Parameters 



Parameter Value 



Mp [Mjup] 0.36 ± 0.03 

Rp [i?jup] 0.71 ± 0.05 

\oggp [cgs] 3.357in3o 

Pp[gcm-3] 1.25 ±0.28 

Semimajor axis, a [AU] 0.0432 ± 0.0006 

e cos 00 -0.0014 ±0.0012 



Method 



Spectroscopic orbit^ 

Rp/Ri, from light curves and i^y^ = 1.45 ± 0.10 
Light curve and spectroscopic orbit^ 
Mp, Rp given above 
Kepler's Law'^ 

Timing of secondary eclipse'^ 



^Using K = 43.3 ± 1.2 m s'^, from Sato et al. (2005); P from Table 2; and = 1.30 ±0.06, 
based on the Yonsei-Yale isochrone fit to a/i?* and T^s. 

^Using K = 43.3 ± 1.2 m s~^, from Sato et al. (2005); P from Table 2; and i, a/Rp from 
the light curve analysis. This method is described in detail by Southworth et al. (2007) and 
Sozzetti et al. (2007). 

'^Using P from Table 2; and = 1.30±0.06, based on the Yonsei-Yale isochrone fit to a/R^ 
and Teff. 

^Using the secondary eclipse time HJD = 2453606.960 ± 0.001 from Harrington et al. (2007) 
and the ephemeris given in Table 2 (after correcting for the 43 s light travel time across the 
orbit). 



